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Abstract. - We investigate the behavior of a two dimensional array of Bose-Einstein condensate 
tubes described by means of a Bose-Hubbard Hamiltonian. Using a Wannier function expansion 
for the wavefunction in each tube, we compute the Bose-Hubbard parameters related to two 
different longitudinal potentials, periodic and quasiperiodic. We predict that - upon increasing 
the external potential strength along the direction of the tubes - the condensate can experience a 
reentrant transition between a Mott insulating phase and the superfluid one. 



Introduction. — Since Greiner et al. first succeeded 
at storing a Bose-Einstein condensate in a two-dimensional 
array of narrow tubes using an optical lattice pQ , experi- 
mental progress on low-dimensional Bose systems has been 
tremendous. To name but a few highlights, the Mott tran- 
sition was observed in three [2J, one [3J, and two dimen- 
sions [HG3, and in ID a Tonks- Girardeau gas has been 
realized 6,7. Moreover, in 2D the Kosterlitz-Thouless 
transition was observed [8j[9]. The dimensional crossover 
between three, two, and one dimensions for bosons in an 
optical lattice was studied theoretically using Tomonaga- 
Luttinger liquid theory in Refs. Q3I1III] and using Monte 
Carlo simulations in Refs. (T2JE5]. These references stud- 
ied 2D or 3D optical lattices in which atoms can tunnel 
easily along one Cartesian direction but not along the oth- 
ers. The system can then be described as an array of tubes 
of bosons, which may or may not be mutually phase co- 
herent. As the tunneling between tubes is varied, such 
a system will undergo a transition from a 3D superfluid 
(3D SF) to a 2D Mott insulating phase (2D MI), which 
consists of a decoupled array of ID tubes. Similarly, if the 
tunneling probability along all three Cartesian directions 
are made unequal, a 2D SF state can be realized, in which 
there is superfluidity within separate 2D layers but no co- 
herence between them. These transitions are present only 
if the tubes have finite length. 

The theoretical approaches used in Refs. [TUr[T5] are able 
to describe phase fluctuations within the tube-like fila- 
ments of bosons, but they do not capture any nonlinear 
effects due to the possible variation in the width of these 
tubes. Gross-Pitaevskii theory describes such variations 



[14]; however, it does so at the expense of not being able 
to account for phase fluctuations. Keeping these limita- 
tions in mind, we offer in this Letter a description that 
is complementary to those of Refs. [l"0rfl"2] . using Gross- 
Pitaevskii theory in order to understand how nonlinear 
effects may affect the phase transitions in this peculiar 
type of optical lattice. 

Such effects become important if the potential barri- 
ers along the strongly coupled direction are so weak that 
each tube can be considered as a quasi-lD Bose-Einstein 
condensate and the number of bosons is large. It will be 
shown that these nonlinear effects can give rise to a reen- 
trant Mott transition in the array of ID tubes. 

The dilute Bose gas in an external potential V ex t(r) is 
described by the second quantized Hamiltonian 
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The physical setting consists of a two-dimensional opti- 
cal lattice in the x and y plane with period d, creating a 
square array of tubes which develop along the z direction. 
Moreover a weaker optical potential parallel to the tubes 
is added, so that in the present case the external potential 
is given by 

V cxt (x, y, z) = -V x cos — V y cos — — + V z (z) . (2) 

As in Ref. [TS] the many-body field operator of the gas 
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V>(r) is rewritten as a sum of local operators 

jxjyk 



(3) 



where (j) ]xjyk (r) = fe (a; + j x d, y + j y d, z) and each b jxjyk 
acts on the k th state of the tube in the position {j x ,j y ). 
The 4>k are a complete set of wavefunctions. In the present 
setting we expect only the lowest energy state to be occu- 
pied, so that it is possible to drop the k index from ([3]). 
Using the expression ([3]) the second quantized grand po- 
tential 



can be rewritten 
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where j — (jx,jy), hj = blbj, and the inter-tube tunneling 
matrix elements are 
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with a — {x, y}. The in-tube interaction energy is 

U = g fdr\0(r)\\ (7) 



jj, = — dr 



— |V0(r)| 2 + |^(r)| 2 y cxt (r) 
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where /zq is a constant that fixes the number N tot of par- 
ticles in the whole system and g = Airh a s /m. 

Finding the wavefunction (f>(r) is a nontrivial task: 
We are ultimately concerned with approximating the full 
many-body ground state of the 3D system. In the ab- 
sence of nonlinear on-site effects, the maximally localized 
Wannier functions can be used as a single-particle basis in 
which to expand the full many-body problem [15JI16] . For 
the problem at hand, we are not only concerned with non- 
linear effects from on-site repulsion, but also the elongated 
nature of the tubes. A rigorous recipe for how to choose 
the basis here cannot easily be formulated, but in prin- 
ciple the accuracy of the chosen basis could be assessed 
a posteriori by estimating how well the final solution ap- 
proximates the 3D many-body ground state. 

Nevertheless, we take guidance from the ID single- 
particle problem and apply the following ansatz to de- 
scribe the wavefunction in each tube, 



( r ) =J2 a ij( z ) w i(x)w 3 (y), 



(9) 
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Fig. 1: Left panel: behavior oit/U as a function of V z for some 
values of g n, using the the periodic potential (|12p with V = 
77.35-Er, as computed with the Wannier function approach. 
Upper right panel: behavior of t/U as a function of V z for gn = 
11.08 Er. Lower right panel: tentative MF order parameter 
(b) against V z (see text) for the case L = 28 (11.9/mi), and 28 
particles for each tube. 



(6) where <fi is normalized to unity 



dr<f>(r) = 1 



(10) 



and Wi(x) and Wj(y) are the Wannier functions for our lat- 
tice configuration along the x and y direction respectively. 
Since we expect the ground state in each tube to be sym- 
metric, the sum employs only the even Wannier functions. 
The energy functional for the Bose-Einstein condensate in 
each tube is 



H[4>,4>*] = dv — |V0(r)| 2 (11) 

J vim 

+ U cxt (r)|0(r)| 2 
+ — Wr) 

with n the mean number of particles in each tube. Insert- 
ing (|9]) in Eq. (jlll) , an expression of the total energy as a 
function of the coefficients a,ij(z) is obtained. The ground 
state has been found by minimizing this expression using 
the conjugate gradient method; in the chosen set of pa- 
rameters it is enough to employ < i,j < 4, since the 
occupation number of higher Wannier states is negligible. 
Once the ground state is obtained, the coefficients t and 
U can be computed by means of ([6} and (jTj) . 

Periodic potential. — Let us suppose that the po- 
tential along the tubes is periodic. For definitcness, let it 
have the same wavelength of the potentials in the x and y 
directions 

V z {z) = V z cos(^z) . (12) 
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Focusing on the case in which V = 77.35 Sr (Er is the 
recoil energy Ti 2 /md), in Fig. [T] (left panel) we plot the 
behavior of t/U while varying V z from to 60 -Er, for 
some values of the interaction strength. We see that for 
gn — 11.08 E R the ratio t/U decreases until V z ~ 39 E R , 
then it begins to rise. It is known [T71HB] that the Mott 
insulating phase appears as a set of "lobes" in the plane 
n/U-t/U, each one corresponding to a precise number of 
particles Atubc for each tube; outside of the 2D MI zone 
the system moves along lines of constant A^ube- The initial 
decrease and subsequent decrease of t/U therefore means 
that for an appropriate length of the tubes the system 
can cross the 3D SF - 2D MI transition twice. In order 
to give an estimate of the critical potential strength, the 
Gutzwiller approximation is employed (see for example 
P23E5]). The condensate order parameter is given by the 
expectation value of the destruction operator (b). FigurefT] 
(right lower panel) plots the behavior of this parameter 
for L = 28 (11.9/zm) and 28 particles in each tube. The 
bounded interval in which the superfluid phase vanishes 
is a clear consequence of the nonlinear behavior of the 
system. 

However, this example of a reentrant phenomenon does 
not survive a more accurate analysis. In fact, the param- 
eters are such that when the system enters the 2D MI 
phase - thus isolating each tube from one another - the 
gas of bosons is in an insulating phase within the single 
tube. The assumption of coherence along the z axis is not 
valid anymore, and a reentrance to a 3D SF is question- 
able. Even changing the parameters, we were not able to 
find a case where the system remains superfluid along the 
tubes in the critical zone. This situation, added to the 
experimental difficulty to obtain the precise combination 
of parameters that realize the reentrant transition, makes 
the described phenomenon very unlikely to be observed in 
the case of a periodic potential along the tubes. 

Quasiperiodic potential. — In order to avoid the 
Mott insulating phase within the tubes we employ instead 
a quasiperiodic potential 



V z (z) = V z cos 
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z ) + V z cos 
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where A = 89/55 is an approximation to the golden ratio. 
We have chosen to use the same lattice parameters as in 
the previous analysis, with the length of the tubes chosen 
to be L = 89 d (37.825/mi). Proceeding in this way we see 
a behavior similar to the one found for the periodic case: 
t/U as a function of V z is a non-monotonic function for 
g n — 8.99 En, as shown in Fig. [3] This situation gives rise 
to a reentrant transition when the number of particles in 
each tube is between 2500 and 2700. Figure [2] (lower right 
panel) shows the mean field order parameter for iVtube = 
2600. According to the Gutzwiller approximation, in the 
shown parameter range there are a maximum of (b) 2 ~ 100 
particles in the condensate. 

In order to better understand the reason behind this 
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Fig. 2: Left panel: behavior of t/U as a function of V z for 
the quasiperiodic case; the picture refers to some values of gn, 
while V — 100.0-Er, as computed with the Wannier function 
approach. Upper right panel: behavior of t/U as a function of 
V z for gn = 8.99 Er. Lower right panel: MF order parame- 
ter (b) against V z for the case L = 89 d (37.825/im) and 2600 
particles per tube. 



phenomenon, the plot in Fig.|3]shows the behavior oft and 
U against V z separately for the case gn = 8.99 Er, and 
Fig. @] shows the density of the wavefunction J dy \<fi\ 2 . In- 
creasing the longitudinal quasiperiodic lattice we see that 
the peaks in the boson gas profile become narrower, thus 
raising the value of U (J7J), while at the same time the 
wavefunction widens in the radial direction, increasing the 
tunneling rate t (JSJ). At first U rises faster than t but for 
t « 25 Eii this relation is reversed. The observed nonlin- 
ear effect is therefore the result of a competition between 
U and t, i.e., between the increased density in a lattice 
minimum and the broadening of the tubes upon raising 
the longitudinal lattice strength. 

In the presence of the quasiperiodic potential that we 
employ, the Mott insulator phase within a single tube is 
inhibited, but a new quantum phase comes along: the Bose 
glass. Recently a Bogoliubov analysis has been employed 
to trace the boundaries of the Bose glass phase |20|21) , and 
this method is expected to be accurate within the limit of 
the Bogoliubov approximation, i.e., weak interaction and 
high density {g/n <C H 2 /m [22]), as is the case for our 
current work. Using the same method we were able to 
determine that the single tube never leaves the superfluid 
phase during the reentrant phenomenon, thus making it 
possible for the system to experience twice the 3D SF- 
2D MI transition while increasing the quasiperiodic lattice 
strength. 

We thus conclude that our model is appropriate to pre- 
dict the phase transition, and the consequent reentrant 
behavior. 

Conclusions. — In this Letter we have studied the be- 
havior of a 2D array of strongly elongated tubes of bosons 
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Fig. 3: Behavior of t and U as a function of V z for gn — 
8.99 £_r, in the quasiperiodic case. While they both grow upon 
increasing Vz, the different rate at which they increase is re- 
sponsible for the nonlinear behavior of t/U (Fig. [2]) 



Fig. 4: Surface plot of the density profile J dy \<j>{x, y, z)\ 2 for 
- from left to right - V z = 2, 10, 25 and 40 £r, shown for < 
z < 10 d. Increasing the longitudinal quasiperiodic lattice we 
see that the peaks in the boson gas profile become narrower, 
and at the same time the width of the wavefunction increases. 



in an optical lattice, employing the Bose-Hubbard Hamil- 
tonian in order to describe the 3D SF - 2D MI transi- 
tion. Using a Wannier function ansatz and minimizing 
the Hamiltonian of the system we were able to compute 
the ground state in each tube, and link the parameter t 
and U of the Bose-Hubbard Hamiltonian to the physical 
parameters of the atom gas. Using this approach we were 
able to see a nonlinear behavior of the parameters, leading 
to the possibility of a reentrant transition between the 3D 
SF and the 2D MI phase. 

We applied our method to two different potentials along 
the tubes: a periodic and a quasiperiodic modulation. 
While the non-monotonic behavior of t/U is seen - for 
a specific choice of the parameters - in both cases, the 
periodic potential puts the tubes in a ID MI phase, 
which would nullify the phenomenon we want to ob- 
serve. In the quasiperiodic case, we predict that in an 
array of 87 Rb tubes 37.825 fim long, with d — 425 nm, 
Vx = V y — 100 Er, interaction strength gn ~ 9 En, and 
Ntubc ~ 2600, the system goes through the 3D SF - 2D MI 
transition twice while increasing V z . The system itself is 
found to be in the insulating phase for V z between ~ 10 Er 
and ~ 30 -Er. We stress that this is not a consequence of 
a localization phenomenon, but the quasiperiodic poten- 
tial is only used to put the tubes in the right parameter 
regime. 

Future prospects include investigating how the algebraic 
decay of phase coherence - which is not included in this 
study - affects the detailed values of the critical parame- 
ters. Moreover, a real experiment will employ an inhomo- 
geneous trapping potential. This introduces the complica- 
tion that not all the tubes would be equally long and the 
occupation number fluctuates. What is more likely to be 
observed in this situation is a dip in the phase coherence 
while varying V z in the range described above. Finally, we 



expect the combined quasiperiodic-plus-periodic potential 
under consideration to have a rich phase diagram if a wider 
parameter regime is considered; this phase diagram will 
include MI in different dimensions, Bose glass, and com- 
binations thereof. The investigation of these phases will 
be the subject of a separate study. 
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